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Abstract. The spectrum of guided modes of linear chain of dielectric and left-handed 
cylinders is analyzed. The structure of eigenfrequences is much more richer if cylinders 
are made from the left-handed material with both permittivity and permeability 
negative. The number of guided modes is much larger, and their interaction with 
incident electromagnetic wave is much stronger. For some value of the wave vector, no 
guided modes were found. We discuss how these specihc properties of guided modes 
correspond to folded bands, observed recently in photonic structures with left-handed 
components. 
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1. Introduction 

Spatial periodicity of photonic structures composed from two dielectric materials is 
responsible for a broad variety of interesting physical phenomena. Band structure of 
frequency spectrum, appearance of gaps of forbidden frequencies EEIE], bound states 
[1], surface states [5] and guided modes in periodic structures P [7| represent only a 
few of them. Artihcial metal - composite with negative effective permittivity can be 
constructed when one material is supplied by metal PPCQ]. 

Discovery of left-handed materials HB m with simultaneously negative 
permittivity and permeability, naturally addresses the question whether such materials 
could be used in the construction of new photonic structures and how negativity’s 
of permittivity and permeability influence properties of resulting composite. Already 
in one dimensional systems, a new gap associated with zero mean value of refractive 
index was observed [13]. In two dimensional photonic crystals, the use of left-handed 
components led to the existence of folded bands with zero density of states for a 
certain interval of wave vectors mi US]- Numerical calculation of the transmission 
coefficients for photonic crystals with left-handed cylinders led to unexpected numerical 
instabilities when standard numerical methods were used [T6|- Typical example of such 
instability is shown below in figure which shows numerical data for transmission 
coefficient obtained by the transfer matrix algorithm for different spatial discretization 
give inconsistent results, which do not converge even in the limit of the most accurate 
accessible numerical approximation. Physically, such numerical instability indicates that 
the spatial distribution of electromagnetic field inside the photonic structure is strongly 
inhomogeneous [HI- 

In this paper we suggest that unusual properties of periodic photonic structures 
with left-handed components are due to excitation of broad spectra of guided modes 
which interact with incident electromagnetic wave. We show that indeed photonic 
structures with left-handed component support the excitation of a broad variety of 
spatially inhomogeneous guided modes which strongly influence the transmission of 
electromagnetic waves. 

The paper is organized as follows: in Section we describe the model and derive 
mathematical equations for guided modes in the most simple photonic structure: the 
linear chain of cylinders. Obtained system of linear equations is used in Section to 
calculate of complete spectrum of guided modes both inside and outside the light cone. 
In Sectionj^we use our numerical method for the calculation of the transmission through 
of plane electromagnetic wave through our structures. For dielectric cylinders, we proved 
that Fano resonance observed in the transmission spectra are due to the excitation of 
leaky guided modes miniiii. Much stronger interaction of transmitted wave with guided 
modes is found in structures containing left-handed cylinders. Conclusion is given in 
Section IH 
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Figure 1. Structure consists from periodic linear chain of homogeneous cylinders 
of infinite length in the z direction. Cylinders have relative permittivity e and 
permeability p. We consider £ = 12 and /j, = 1 for the dielectric cylinder, and £ = —12 
and /i = — 1 when cylinders are made from the left-handed material. The chain is 
located along the line ?/ = 0. Central cylinder is located a.i x,y = 0,0. The structure 
is periodic in the x direction with period a. The radius of cylinders is i? = 0.3a. 
Electromagnetic wave impinges on the structure in y direction with polarization E || z 
or iJ 11 z. Since the structure is periodic in the x direction, it is sufficient to consider 
only one unit cell. 


2. The model 

The structure of interest consists of an infinite periodic chain of cylinders embedded 
in the vacuum, (hgure Distance between cylinders, a, dehnes spatial periodicity of 
the structure along the x direction. Cylinders are inhnite along the z; direction, have 
radius R = 0.3a and are made from homogeneous material with relative permittivity e 
and permeability /i. We chose £ = +12, /i = +1 for dielectric cylinders, and £ = —12, 
yU = — 1 for cylinders made from left-handed material. 

We are interested in eigenstates propagating along the x direction with wave vector 
q and frequency u. If {q,uj) lies inside the light cone (leaky modes with q < u/c), these 
states could couple to incident electromagnetic field. 

2.1. Equations for guided modes 

Consider hrst the polarized electromagnetic wave. For the central cylinder, we 
express the electric intensity and tangential component of the magnetic intensity ht 
in terms of Bessel functions [TH]. Inside the cylinder (r < R) we hnd 



fc >0 


fc >0 


( 1 ) 


iZEf{r,ip) = + 2^afXcos{kip) + ‘2i'^a^Xsm{kip) 


k>0 


k>0 


In equation [T] we used 


Jk = Jki^TTnr/X), X = J'^{2'Knr / X) 


( 2 ) 
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where Jk, J'k are the Bessel function (fc = 0,1,...) and their derivatives, A = ‘I'Koju is 
the wavelength of electromagnetic field in vacuum, n = ^/ejl is the index of refraction of 
the cylinder and Z = pi/e (c is the light velocity in vacuum and Zq is the vacuum 
impedance). Unknown parameters a'^ and a~ determines amplitudes of even (oc cos k^p) 
and odd (oc sin kip) cylindrical waves, respectively. 

The electric and magnetic fields outside the central cylinder consist from two 
contributions: the first one is the field scattered on the central cylinder, 

e°(r, (p) = 'Ho/Sc)’ + cos(fc(p) + 2iy^p^Hk sin(fc(p) 

fc >0 fc >0 

( 3 ) 

-iZoA?(r,V>) = WM + 'i'Z^K'S'kCoaikip) + 

k>0 k>0 

(r > R) with free parameters f3~^ and /3“. Here, 

Rk = Hk{2nr/X), U/= H/{2Tir / X). ( 4 ) 

and Hk{z) = Jk{z) + iYk{z) is the first Hankel function [THl fTU] . 

The second contribution to the external fields is represented by fields scattered 
from all other cylinders. For the nth one (n = ±1, ±2,..., Ng) we express electric and 
magnetic intensity in coordinates associated with the center of the cylinder (figure 

= Ho{Wn)l3fo + 2^Pfi.Hk{Wn)cOs{ken) 

k>0 

+ /3^f^Hk{wn) sm{k9n) 

k>0 


-i^ohZ(^n,9n) = ffo(uJn)/3fo + 2^l3f^^H/{wn)cos{kdn) 

k>0 

+ ‘^^'^^nkHkiWn) sm{k9n) 


( 5 ) 


Zoh/Jfn,9n) = Ho{Wn)/3fQ + 2 ^/ 3 +^— Hk{Wn) cos{k9n) 


k>0 


_ 1 ^ 

- HkiWn) sm{k9n) 


k>0 

where we use 

Wn = 2 'Kf,n/X ( 6 ) 

Thanks to the periodicity of the model, coefficients of the nth cylinder could 
be expressed in terms of coefficients as 

^nk = ( 7 ) 

Coefficients a and f3 could be calculated from the requirement of the continuity of 
fields Cz and h^ at the surface of the cylinder: 


( 8 ) 
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Figure 2. Parameters used in derivation of equation 


and 


hf{R ) = h'f^^ = On) cosa„ - On) sin(9) 


n^O 


where an = On — ‘P- 

Before solving this system of equations, we hrst transform all helds in equation 
as a functions of R and p. This can be done with the use of the Gegenbauer formula 
for cylindrical functions (hgure 

+ 00 

Hn,{w)e^^^^ = Y, Hn,+k{u)Jk{v)e^^^^ ( 10 ) 

k = — OQ 

[19] and similar formula for the hrst derivative H'n^{w)\ 


+ 00 


( 11 ) 


k=—oo 

Inserting ( [II| ) into equations ([^ we express, after some algebra, the helds at the outer 
boundary of the cylinder in the form 

N N 


e°f^\R+,p) = Y^km/d:^coskp 

k,m=0 

N 


Y CkruPm sin kp 

'T" ( 12 ) 


h°"*(i?+,(p) = cosfcip + Y. C\nildm^W-kp 

k,m=0 k,m=l 

Where we introduced N as the highest order Bessel function used in the expressions of 
the helds. Explicit form of matrices B, B’, C and C’ is given in the Appendix A. 



Figure 3. Parameters used in Gegenbauer’s relation, equations [T0| and [m 
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Incident electromagnetic field contributes to the rhs of equation We express it 
in the cylindrical coordinates 

N N 

= Mv) + cos(/c(p) + i'^e^ sm{k^) 

k>0 k>0 {-\ n\ 

N N 

-iZoh^{r,^) = J[^{v) + ^hfcos{kip) + i^h^ sm{kip) 

k>0 k>0 

Inserting fields ([^ and 12) into equations and we obtain the set of linear 
equations for even modes: 

+ ^k 

(14) 

CAW = E™B’,^/3+ + h+ 

(fc = 0,1,..., N) and similar equations for odd modes 
^kOCk ~ ^kmPm + Cfc 

(15) 

CA kOif. = 'Yhm ^ km/djn + ^k 

{k = 1,2,..., N). Here, C = Z/Zq, Ak = 54(2 - 5ko) and A4 = Jl{2 - 6ko) (equation 

0 . 


In the next step, we remove a'^ from equations [1^ and obtain the system of linear 
equations for parameters 


E 


^km S» rrf^ km 
^k 


/3+ = e+- C^ht, 

'-^k 


and equivalent equation for parameters fd : 

E 


\Jk 

_ 

km Sj rr( km 
*^k 


Pm = Cfc -C^h^, 

’-tk 


k = 11,1,..., N 


k = l,2,...,N 


(16) 


(17) 


With known parameters fd^, we find from equations |14| and |16 


± 


In numerical analysis, we consider the number of cylinders W = 100 — 4000 and 
number of modes N < 22. We found that consideration of iV = 12 modes is sufficient 
in most of analyzed problems. 


2.2. Transmission coefficient 

Transmission coefficient can be calculated as the ratio of the ^/-component of the 
Poynting vector Sx{yp), calculated for any Pp > r to the incident Poynting vector, 

SI 


T = 


SxiVp) 

Si 


r-\-aj2 r-\-aj2 

/ dxe,{x,yp)hl{x,yp) / / dx el{x,yp){hl{x,yp))*{18) 
'-a/2 J-a/2 
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Figure 4. Guided modes of the chain of dielectric cylinders {e = 12, ^ = 1). Left 
figure shows the frequency dependence of diagonal elements In A^k for <7 = 0 (equations 
and Four even modes and three odd modes can be identified. Right figures 
presents the spectrum of Ez polarized Open and closed symbols correspond to even 
modes (n = 0 and 1) and odd mode n = 1. Owing to the periodicity of the problem, 
it is sufficient to consider wave vectors 0 < q < ir/a [T]. 


2.3. Polarization 


The above analysis Can be repeated for the polarized modes. It turns out that hnal 
equations (16) and 0 remain the same if we substitute 




(19) 


and exchange of Cx and hz in expression for the transmission coefficient (18). 


3. Guided modes 


Consider first equations IT and 17_ without incident electromagnetic wave. Then, the 
spectrum of guided modes can be found from the requirements of zero determinant 


and 


det 


det 


■p _ 

^km S 

^k 


■'km 


3^k 


= 0 


= 0 


for even and odd modes, respectively. 

With the use of Gauss-Jordan elimination method 


( 20 ) 


( 21 ) 


we transform matrices in 


equations (20) and (21) to diagonal form and plot the frequency dependence of inverse of 


obtained diagonal elements A^k- This enables us not only to hnd the eigenfrequency of 
guided modes and their lifetime ED. but also identify their symmetry. As an example, 
we show in the left figure 4 resonances of leaky guided modes with zero wave vector 
g = 0. 
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Figure 5. Spectrum of guided modes of the chain of cylinders made from the 
left-handed material (e = —12, /i = — 1) cylinders. There are no modes around 
Qc « 0.24 X 27r/a. Open and closed symbols correspond to even and odd modes, 
respectively. Note that guided modes changes symmetry when q crosses q^. 


Right figure 4 presents eigenfrequences of guided modes as a function of wave vector 
q identified from the position of resonances in corresponding diagonal elements. We 
found two even modes with k = 0 and 1 and one odd mode k = 1, which agree with our 
expectation [1]. 

The spectrum of guided modes for left-handed cylinders is more complicated and 
consists of series of even and odd modes. In contrast to dielectric cylinders, the 
eigenfrequences of guided modes decreases when mode index k increases. The number 
of modes depends on the model parameters and increases when absolute value of the 
refractive index increases. 

Another unexpected property of guided modes is the existence of ’’critical value” 
of the wave vector [q^ ~ 0.247 x 27r/a in our model) for which no guided mode exists 
(figure 5). Preliminary analysis of other left-handed structures with p = —1 (data not 
shown) indicates that q^ depends neither on (negative) permittivity nor on the radius 
of cylinder. 


4. Interaction of external field with guided modes 


We solve equations IT and 17_ and calculate transmission coefficient for incident plane 
wave with wavelength A propagating along the y direction. First, we express the 


± 


coefficients and 


~ ~ + (~1)^)/(1 +'^fco) k — 0,1,... 


( 22 ) 


et = 'h = All'll - (-1) ). 4=1,2,... 

{v = 2 'kR/\). From known coefficients (3^ and /d“ we find spatial distribution of fields 
and calculate transmission coefficient given by equation 
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Figure 6. Transmission of the Ez polarized plane wave propagating through an 
infinite ID chain of dielectric cylinders. Top panel displays the transmission coefficient 
as a function of a/A, calculated by the transfer matrix method and by the formula 
(18). Three Fano resonances correspond to maxims of coefficients /3 shown in middle 
panel. Real part of the intensity of electric field inside the cylinder for three resonance 
frequencies, a/A = 0.549, 0.759 and 0.958 is shown in the bottom panel. 


4-.1. Dielectric cylinders 

Figure shows the transmission coefficient for the polarized plane electromagnetic 
wave as a function of wavelength. Obtained results agree very well with numerical 
data calculated by the transfer matrix method [23l [23] . Typical Fano resonances in the 
transmission spectra lac] are clearly visible and could be associated with excitation 
of corresponding guided modes [22] in cylinder chain. Owing to the symmetry of the 
incident wave which does not depend on x, only modes symmetric with respect to 
X —>• —X could be excited. Two even modes with k = 2 and fc = 4 and one odd mode 
{k = 3) have been excited. Bottom panels show the real part of electric held inside the 
cylinder and con&m predicted symmetry of excited states. 

Similar results have been obtained also for the Hz polarized incident wave shown 
in hgure 7. In this case, two resonances of A; = 0 mode are clearly visible. The lower 
one, together with the hrst odd resonance, deform strongly the transmission spectrum. 
Three resonances observed at higher frequencies could be easily detected both from the 
transmission spectra and from the frequency dependence of parameters (d. 
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Figure 7. Transmission of the polarized plane wave propagating through an 
infinite chain of dielectric cylinders. Top panel displays the transmission coefficient 
as a function of a/A, calculated by the transfer matrix method and by the formula 
equation 1 18| Fano resonances correspond to maxims of coefficients j3 shown in middle 
panel. Real part of the intensity of magnetic field inside the cylinder for three resonance 
frequencies, a/A = 0.729, 0.865 and 0.931 is shown in the bottom panel. 


4 . 2 . Left-handed cylinders 

The same analysis of the transmission spectra of cylinders made from the left-handed 
medium (e: = —12, /i = —1) is more difficult, since standard numerical techniques 
(transfer matrix method, RCWA ) seem not to be suitable for the calculation of the 
transmission coefficient. As shown in the left panel of figure results obtained by 
the transfer matrix method with three discretization of the unit cell provide us with 
completely different transmission coefficients. Such failure, obtained also when the 
transmission was calculated by the RCWA method indicates that the eigenstates of the 
electromagnetic field inside the structure are strongly inhomogeneous. Consequently, 
excited guided modes cannot be expanded in a finite series in the the basis of 
eigenfrequences, used in the transfer matrix method. 

The assumption of strong spatial inhomogeneity is confirmed by frequency 
dependence of coefficients ft shown in middle panel of figure Incident electromagnetic 
wave excites in the left-handed chain a series of guided modes. The resonant frequencies 
of these modes lie very close to each other (in fact, since the resonances have a finite 
width, they overlap). Also, absolute values of coefficients ft are in order of magnitude 
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Figure 8. Left: transmission of the polarized electromagnetic wave through the 
linear chain of left-handed cylinders calculated by the transfer-matrix method with 
three different space discretization (120, 240 and 360 mesh points per unit cell). Middle: 
spectrum of coefficients /3. Note the scale of the vertical axis. Right panel shows the 
transmission coefficient calculated by the present method (equation 18 ). Bottom 
panel shows real part of the intensity of electric field inside the left-handed cylinder for 
frequencies associated with resonances of even modes (a/A = 0.027, 0.078 and 0.184), 
and odd modes (0.016, 0.047 and 0.140). 


larger than in the case of dielectric cylinder. 

Right panel of figure]^ shows the transmission coefficient obtained by eqnation (18). 
Two minims in the transmission coefficient coincide with eigenfreqnences of excitation 
of gnided modes. The spatial distribution of electric field shown in the bottom panel 
conhrm that many resonances are excited simultaneously - in contrast to dielectric 
structure, we cannot easily identify the order of excited modes. 


The above mentioned problem of numerical instability is not actual when 
transmission of the polarized wave is calculated (figure 9). Now, numerical data 
obtained by the transfer matrix method agree perfectly with those found by the present 
method (left panel). As shown in the right panel, resonances of guided modes are very 
weak and broad, so that they have negligible influence on the transmission coefficient. 
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Figure 9. Left: Transmission of iL^-polarized plane wave through the chain of left- 
handed cylinders calculated by the transfer matrix method and from equation |18| 
Similarly to the case of dielectric cylinders, two numerical methods give identical result, 
because Fano resonances (shown in the right panel) are much weaker than in the case 
of Ez polarized wave (Fig. |^. 


5. Conclusion 

We presented physical and numerical analysis of the guided modes of linear chain of 
cylinders, made either from the dielectric or from the left-handed material. Comparison 
of these two structures shows new phenomena could be observed when permittivity 
and permeability possess negative values. We proved that the transmission of 
electromagnetic waves through the photonic structure with left-handed components is 
strongly influenced by guided modes excited in the structure. 

The spectrum of guided modes contains much more branches. In contrast to 
dielectric structure, the eigenfrequency of these modes decreases when mode index k 
increases. Also, there are no guided modes for the wave vector close to some critical 
value Qc- We suppose that this absence of guided modes corresponds to the folding of 
frequency bands observed recently HSIE]. 

Very rich spectrum of guided modes is responsible for numerical instabilities when 
the transmission coefficients is calculated by standard methods. Since eigenmodes of 
the left-handed structure strongly differs from plane waves, any numerical algorithm, 
based on the expansion of the fields into the plane waves (transfer matrix or RCWA) 
fails to recover true transmission coefficient. 
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For k,m = 1,2,... Nb, Jk = and Rk given by eqnation we obtain the 

A + 1 X A + 1 matrix B 


Ns 


-Boo 

= Ro¬ 

f ^ 2 cos 

qHo{un)Jo 


ns 

n=l 


Bom 

= V 

Hmiun)Jo 

X [(_X)”^g*'?“’^ -|- Q-iqan- 


n=l 

Ns 



Bko 

= v 

n=l 

Hk{un)Jk 

X 2[6*9“^ + (-1)^6-*''“”; 


(23) 


Ns 


Bkm = 2nk6km + J2 + + 

n=l 

and the A x A matrix C; 

Ns 

Ckm = ^RkSkm + 


n=l 


(24) 


Matrices B’ and C’ conld be obtained from B and C, respectively, by snbstitntions 

Jk — t J'k, R-k — t R'k (25) 
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